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A model for precipitation of the plate-shaped second-phase under applied stress is presented. The 
precipitates in the matrix-precipitate system are represented by their local volume fraction ^ and 
an orientation parameter 6 that defines the alignment of a precipitate platelet in a given direction. 
Kinetic equations, based on diffusion theory and classical nucleation theory, are used to describe 
the time evolution of and 9 . The model is used to describe the stress orientation of hydrides in 
Zr-alloys in light of experiments. 

PACS numbers: 



Stress-induced orientation of second-phase particles occurring during precipitation is ubiquitous in many metallur- 
gical systems. Stress orientation is referred to the selective precipitation of a second-phase on a preferred habit plane, 
caused by the application of an external stress during aging or cooling. For example, in the Ti-H system, hydride 
precipitation has been observed to occur under an applied tensile stress, upon w^hich titanium hydrides form on the 
titanium habit plane oriented most nearly perpendicular to the stress axis Similarly, in polycrystalline zirconium 
alloys, it has been observed that stress-orientation of hydrides occurs under tensile stresses, where hydride platelets 
would get aligned perpendicular to the stress axis and can lie in specific habit-planes; namely on either {lOlZ} or 
{112Z} crystal planes 0. Under compressive stresses, however, hydride platelets tend to orient parallel to the direction 
of apphed stress Q. 

The stress orientation effect has also been observed in Fe-N alloys, in which FejfiNo platelets {a" phase) precipitate 
on a specific habit plane perpendicular to the direction of applied tensile stress d, [5| . Other alloy systems exhibiting 
stress orientation effects include Ni-base alloys with Ni-Nb precipitates [3| , Fe-Mo- Au with Au particles [3] , and Al-Cu 
alloys with 9' second-phase particles ^ ^, JJ)J . Similarly, favoured orientation of precipitates has been observed in 
ferromagnetic materials, when the material is cooled or heat-treated below its Curie temperature in magnetic field. 
For instance, magnetic field orienting of Fei6N2 platelets in Fe-N solid solution and alignment of disc- shap ed nitrogen 
atom clusters on the {001} matrix planes of a nitrogen ferrite alloy in the presence of magnetic field [ill fl^ [T^. 




The precipitation of titanium hydride in titanium and zirconium hydride in zirconium alloys indicates that the stress 
orientation of the hydride platelets occurs through orientation of nuclei rather than through selective growth after the 
critical nucleus size has been reached [H, H, H, [3 El) [III • Likewise, quantitative metallographic studies indicate that 
^'-phase alignment in Al-Cu alloys occurs during nucleation and not under the subsequent growth [^, [loj . 

In this paper, we present a model for kinetics of precipitation of second phase in solid solution under an applied 
stress. The model may be considered as a mean field description of the matrix-precipitate system, in which the 
plate-like precipitates are characterized by their volume fraction and a parameter describing their orientation in the 
matrix. We further assume that second-phase orientation occurs during the nucleation stage of the precipitation 
process. Section [ll] describes the model for precipitation kinetics and stress orientation. In section IIIIl a case of 
hydride precipitation and orientation in a zirconium alloy is considered, where the material-dependent parameters 
entering the model are quantified in light of experimental data. The paper concludes with a discussion of results. 



A simple geometrical model for the plate-shaped precipitates in solid solution is envisioned. The precipitates 
consist of a regular array of parallel ellipsoids embedded in a metal matrix. The sample material, initially at some 



I. INTRODUCTION 




II. MODEL 



A. Precipitation kinetics 



* A portion of this paper was presented in MS&T'08, October 5-9, 2008, Pittsburgh, Pennsylvania, USA. 
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high temperature, is cooled below the solubility line. Precipitation proceeds by diffusion of solute atoms from the 
metal matrix to locations, e.g. dislocation lines and/or grain boundaries, where second-phase particles form. A model 
analogous to that used by Ham for describing precipitation of solute atoms on dislocations is considered (see 
appendix E]). In this model, the precipitated fraction of excess solute w can be calculated as a function of time t fairly 
accurately by the simple relationship of the form: w = 1 — e"*/"^", where tq — ^^/(apD), i is the inter-prccipitate 
distance, ao the lowest eigenvalue of the prevailing diffusion equation, and D is the solute diffusivity. Here w is a 
solution of the kinetic equation 

^ + -(--i) = o. (1) 

at To 

It is worth noting that Eg . pertains to a class of kinetic equations used to describe the overall precipitation 
process of excess solute [18|. Now in our convention, we put w — C/Ce, where Ce is the precipitate volume fraction in 
equilibrium, i.e. after infinite time if the current temperature and stress state were held constant. Hence Eq. ^ is 
transformed to 

§ + -(C-C.)-o. (2) 

at To 

which gives the rate of change for the precipitate volume fraction at a given point in the material. Equation ([2]) can 
be solved, subject to an initial condition, (.{t = 0) = Co: and the condition that and tq are constant, yielding 

CW=Ce[l-(l-|)e~*/^«]. (3) 

We consider now a material in which the solubility limit for precipitation is identical to the solubility for dissolution. 
The solute solubility is thus assumed to be independent of the heating/cooling history and depends only on the current 
state of the material. We should note that is a unique function of concentration, temperature T and pressure p in 
the considered material, and can be calculated from the phase diagram by use of the lever rule 

where C is the total solute concentration in the material and Cl and Cjj denote the respective lower and upper 
boundary of the mixed phase region of the phase diagram (temperature vs. concentration). Hence, Ci, is the solubility 
limit in the matrix phase. Equation ([4]) is valid for Cl < C < Cu- It follows that our assumption of constant in 
deriving Eq. ([3]) implies that temperature, pressure and solute total concentration are supposed to be non-varying. 



B. Stress orientation 



We consider a material with precipitate volume fraction and precipitate mean orientation characterized by the 
parameters 0i, which are the fractions of platelets aligned with their principal axis along the coordinate directions 
Xi] see Fig. [T]for a two-dimensional representation. The volume fraction of platelets aligned with the i:th coordinate 
axis, Ci, is then simply Q = C,6i. Hence, its time derivative is 



Next, we make the ansatz 



!^^^^!^+6l— (5) 
dt dt dt 



where rii is the fraction of platelets nucleated in the z-direction at a given instant, rii is a function of the current 
temperature and stress state, and it also depends on the density of precipitate habit planes aligned in the z:th direction. 
The latter quantity is linked to microstructural properties of the material, e.g. the texture and grain shape. Combining 
Eqs. (O and ^ gives 

de^^ru^d^ .^d_^ (7) 
dt C dt dt ^ ' 
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Equation ([7]) provides the evolution laws for the precipitate orientation parameters 9i. It should be emphasized that 
dOi/dt ^ only when dC^/dt > 0, i.e. the precipitate mean orientation is assumed to change as a result of nucleation 
only. If precipitates dissolve, then d^/dt < 0, and 9i are unaffected. This means that the rate of dissolution is 
supposed to be the same for all precipitates, irrespective of their orientation. 

From classical nucleation theory, the fraction of platelets nucleated in the i-direction is calculated as ([l^ and in 
particular appendix to [20]) 

n,(r, aki) = (l + m,o g-'^^^*'^-^^-) (8) 

where rriio = (1 ~ "•io)/"-io; '^lO is the fraction of platelets nucleated in the i-direction under stress-free conditions, 
fl* the critical volume for nucleation of a new phase, /3 = (kBT)~^, T the absolute temperature, ks the Boltzmann 
constant, a^i the stress tensor, and Ae|"; the differential misfit (transformation) strain tensor. For an oblate spheroid 
with major and minor radii a and c, fl* = 47ra*^c*/3, where the asterisk symbol signifies the critical values for 
nucleation. 




0j^l, e^^O dj = 0, 02=1 6'^ = 1/2, 0^ = 1/2 



FIG. 1: Two-dimensional illustration of the variables 9i, which define the mean orientation of platelets. Hence, 9\ = 62 = 1/2 
corresponds to a random platelet orientation in the a::iX2-plane. 

Under constant temperature and stress, Eq. ([7]) can be solved, subject to an initial condition 9i{t = 0) = 9iQ, 
yielding 

where we utilized Eq. ([3]) with r = tq and is 9i at equilibrium, viz. 

Co 

9ie = - {rii - 9io) — . (10) 

In a plane, ^2 — 1 ^ ^ij a-nd therefore a single differential equation is sufhcient to describe the change in the precipitate 
mean orientation. 

Let us now consider a tubular geometry in which the precipitate orientation is referred to the polar coordinates 
(r, if), with 9 = {9r, 9^) and n = (n^, Uip). The differential misfit (transformation) strain tensor is then defined through 

A^i^ell^elr, (11) 

where the right-hand-side terms are the unconstrained transformation strain tensors for radially and circumferentially 
aligned precipitates, respectively. We note that the fraction of precipitates nucleated in the radial direction under 
stress-free conditions, n^o, is related to the fraction of precipitate habit planes parallel to the radial direction. This 
depends on microstructural properties of the material, such as crystallographic texture and grain shape. For Eq. ([8]) 
to be non-trivial, we require that < n^o < 1- 
Let us write Eq. ([8]) in a compact form 



(l + e-*-)"\ (12) 



where $i = v4„(Tfe; Ae^^/T — In (rriio) is a scaled energy parameter. In Fig. [21 we have plotted 9r as a function of the 
reduced time {t/r) and parameter $i using Eqs. ([5|l- (fT^ . In this calculation, we have assumed Co/Ce — 0.025 and 
9,0 = 0. 
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FIG. 2: Generic (scaled) plots of platelet orientation vs. time at different energy scales "1?, defined in the text, using Eqs. 

o-ini). 



TABLE I: Model parameters for Zr2.5Nb alloy, where T is the absolute temperature. 



Parameter unit source 

Ct = 8.080 X 10*exp[-4151.8/r] wppm [21J 
C£ = 2.473 X 10''exp[-3107.9/r] wppm [21J 
Cu = 16641 - 37465 exp[-2414/T] wppm [22] 
Da = 1.17 X 10~^exp[~4041.2/r] m^s [23] 



III. PRECIPITATION OF HYDRIDES IN ZIRCONIUM ALLOYS 



Past theoretical studies on the effect of hydride reorientation in zirconium alloys have shown that the preferential 
orientation of hydride precipitates under stress is most effective during the nucleation stage of the precipitation 
[3, [l^ . It is argued that the driving force for orienting under stress can be a substantial fraction of the overall force 
during nucleation. Furthermore, it is indicated that hydride growth is unlikely to play a role in preferential orientation 

Mm 

Let us first identify material-dependent parameters in the model described in the foregoing section for the case of 
hydrogen in a Zr-alloy, namely Zr-2.5wt% Nb (Zr-2.5Nb), The phases in the binary Zr-H system are as following 
[22| : There are two allotropic forms of zirconium, namely hexagonal closed-pack (hep) a-Zr and body-centered cubic 
(bcc) /3-Zr, two stable hydride phases, face-centered cubic (fee) 5-hydride (« ZrHi.5) and face-centered tetragonal 
(fct) £-hydride (« ZrH2). In addition, there is one metastable hydride phase, 7-hydride, with fct structure (« ZrH), 
which can exist at the lower temperatures in the {a + (5)-phase region. The (S-hydrides appear usually as plate-like 
(oblate spheroid) particles, whereas the 7-hydrides are commonly observed as needle-like ellipsoidal particles. 

We first regard the equilibrium hydride volume fraction as expressed by Eq. Two material-dependent 

parameters appear in this equation, i.e. Cl and Cij. They are considered to be functions of temperature only (table 
m, despite a slight pressure dependence of these quantities. We should also note that Cl differs for hydride dissolution 
Cf and hydride precipitation C£ [2lj. The i5-line phase boundary Cu was obtained by an exponential fitting to the 
data for the H-Zr system presented in We next estimate the relaxation time r in Eq. ([5]). As noted earlier, 

T cx P / Da where ^ is the observed distance between hydrides and Da the hydrogen diffusivity in a-phase zirconium 
(table IJ). For small hydride volume fractions, € can be related to platelet thickness h via [1^ 



_ h(Vrn 

2\ C 



Ih 



(13) 
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where Vm and ( are the volume fractions of the metal and hydride, respectively. Moreover, in a Zr-H alloy, Kearns |26| 
observed that for hydrogen concentrations up to 1000 weight parts per million (wppm), i = 6.7 x 10'^ /i/C, where C is 
wppm hydrogen. With this relation, the relaxation time is determined from measured values of h rather than £. For 
example, for C = 100 wppm, h — 0.5 /im (typical observed value), T = 580 K and Da given in table |T1 we calculate 
r ~ 10 s. Kearns |26| found that a dissolution kinetics model similar to that in Eq. ([2]) fitted well his experimental 
data. He measured the time needed for complete dissolution of hydrides at constant temperatures, ranging from 570 
to 690 K. The hydrogen content of his Zircaloy-4 samples was from 60 to 205 wppm. 

To calculate the fraction of hydrides nucleated in a given direction using Eq. ^ , we need the components of the 
differential misfit strain tensor in Eq. pT|) . Carpenter [27| calculated the unconstrained misfit strains associated 
with the formation of S- and 7-hydride in a-phase zirconium, and expressed the strain components with respect to 
crystallographic orientations of the zirconium matrix; see table |lll Since (5-hydride platelets form on planes, which 
are roughly aligned with the zirconium basal plane, the [0001] direction in table [iTl corresponds approximately to 
the platelet thickness (normal) direction. Hence, the misfit strain is approximately = 0.072 normal to (5-hydride 
platelets, and ej = 0.0458 parallel to the platelets. With these approximations, it follows from Fig. [3] that the misfit 
strain components for radial hydrides are 



Tr _ T ,Tr _ T ,Tr_ T /..N 



Likewise, for circumferential hydrides 



^Tv _ T Tip _ T Tip _ T r\ 



Accordingly, the components of the differential misfit strain tensor in Eq. pT|) become: Ae^^ = ^jf ~ ^± = —0.0262, 

AeJ^ = - ej[ 0.0262, and Ae^^ ^ ej - e"^ ^ 0. Introducing Ae^ = - ej[ 0.0262 and these relations into Eq. 

(|12p . we find rir = (1 + exp[— with $r = il* Aeq {a^^p — arr)/kBT — In(TOro)- This relation for rir can be used 
for cylindrical geometries to calculate the fraction of J-hydride platelets that nucleate with a radial orientation, given 
the stress state and temperature in the material. We may treat m^o and Q* as material dependent input parameters 
to the model, which can be determined from hydride reorientation experiments on a particular material. For example, 
Hardie and Shanahan's (HS) stress reorientation experiment on Zr-2.5Nb containing 100 wppm hydrogen (cf. Fig. 8 
of [1^, samples cooled from 400°C) gives il* /ks = 1.212 x 10^^ KPa~^ and m,.o = 1041. Employing these constants 
for cr^y = 150 MPa, (!„ = 0, T = 573 K, we obtain $ « 1.4 (cf. Fig. We should note that Eq. ^ is formulated 
in terms of true stress, which means that the stress tensor aki comprises contributions from applied loads as well 
as possible residual stresses, thermally induced stresses, etc. Hydride reorientation experiments show that residual 
stresses can have a significant impact on hydride orientation [281 [29j . 



TABLE IL Zirconium hydride misfit strains, as reported by Carpenter [27j . 



Direction in a-phase 
zirconium matrix 


Unconstrained misfit strain 
(5-hydride 7-hydride 


[0001] 


0.0720 


0.0570 


[1120] 


0.0458 


0.0055 


[1100] 


0.0458 


0.0564 



We now rewrite explicitly the expression for the number of platelets nucleated in the radial direction, Eq. ([8]), for 
0, namely 

rir = (1 + mroexp[-(3Q*AeQa^ip]j . (16) 

It is important to point out that the supposition that the parameter fl* is a material constant at best is a rough 
approximation. There exists a supercooling effect, which shows that the formation of a new phase occurs at a 
temperature below the equilibrium solvus temperature Tg by a deviation AT — Te — T. In other words, the critical 
volume SI* or the Gibbs energy of formation G* are decreasing functions of AT — Tg — T (see e.g. [13] )• Furthermore, 
it can be shown that for a plate-like oblate spheroid precipitate (see appendix to [20] ) 

n* = 97r7/7e5;:', (17) 

and the critical energy for nucleation 



G* = ^7/7e^5„-^ (18) 
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FIG. 3: Radial and circumferential hydrides. The hydride transformation (misfit) strain is perpendicular to the hydride 
platelets, and eT parallel to the platelets. 



with the nuclcation energy density expressed as 

9n = ksTCx \n{Cs/Ceq) - 9c + gex- (19) 

Here Cs, Ceq are the concentrations of the solute (hydrogen) in the precipitate (hydride) , in the supersaturated 
matrix (Zr), and in the matrix in equilibrium with the precipitate, respectively. Also, is the misfit strain energy 
(self-energy) per unit volume, and 7/, 7e are the specific interface energies of the flat and the edge side of the platelet, 
respectively |3lj . The interfacial energy may be estimated from measured data on nuclcation rate (3^ . but to our 
knowledge such data are unavailable for the hydride-Zr system. Furthermore, we should note that hydrogen in the 
matrix depletes during the formation of a new phase and therefore Cs is a time-dependent variable. More precisely, we 
may relate Cs to the hydride volume fraction and the total hydrogen concentration in the material C in the manner 

C = Cs{l-C) + CCu. (20) 

where C is governed by Eq. ([2]). Finally, for the calculation of the relaxation time r, we relate the inter-hydride 
distance in the specimen to the volume fraction of hydride by assuming 

(21) 

^i+C(^o-^i)' ^ ^ 

where ^0 and ^l are observational constants. Thus to account for the effect of supercooling, the two ordinary differential 
equations ^ and ([7]) need to be solved simultaneously, by making use of Eqs. @ and (fTB|) - (PT|) . 

Using the aforementioned method, we attempt to simulate an experiment by Hardie and Shanahan [l5| on stress 
orientation of hydrides in Zr-2.5Nb specimens under non-isothermal conditions. More specifically, we consider a case 
in HS's series C tests where the specimen contained C = 100 wppm hydrogen. The specimen was heated from room 
temperature to 673 K to dissolve the hydrides, then cooled under a constant stress (150 MPa) at rates 3.3 Ks~^ 
(673 to 573 K), 2.5 Ks~i (573 to 473 K) and 1.7 Ks^i (473 to 373 K). In our computations, the variable Cl in Eq. 
(ID) is the hydrogen solubility for precipitation Cp. Our analyses show that the choice of Cl has a strong impact on 
the calculated results, since it determines the onset of hydride precipitation during a cooling sequence. Therefore, 
we assume that Cl = SfCp, where 5'/ is a scaling factor, 1 < 5/ < 2, and Cp is given by the expression for C£ 
in table [J We have used the following set of values for the model parameters: gex = 0.0262cr;^;^, = 100 MJm"'^, 
7^ = 0.065 Jm~^, 7e = 0.28 Jm~^, Iq — 100 ^m, li = 0.1 /im, Sf — 1.47, and n^o « 5 x 10~^. In addition, we 
consider (5-hydride with Cx = xNA/Vh, x — 1.66, the molar volume Vh — 1.63 x 10~^ m'^mol^^, and A''^ the Avogadro 
constant. Moreover, we have assumed Ceq = C'[{T) as given in table HI 

We have solved the system of the aforementioned equations by using the Runge-Kutta algorithm of order 4 and 
5 [3^. The results for the time variation of the hydride volume fraction and the orientation parameter 9r are 
depicted in Fig. 21 It is seen that nuclcation of hydride occurs immediately around 3000 s after the start of cooling, 
corresponding to the temperature of about 525 K. Then 9r, after a sharp dip to a shallow minimum, raises and falls 
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slowly to a near equilibrium value of Or = -5, at t = 8000 s, T = 373 K. This value is close to the experimental 
result of Hardie and Shanahan (cf. Fig. 8 of [3). We have also tabulated the results of computations at T = 373 K 
for several initial (total) hydrogen concentrations: C — [100, 200, 300] wppm, and apphed stresses: a^^ — 100 — 300 
MPa. Table [nil shows these results for 9r- The second column in this table can be compared with HS's experimental 
data (Fig. 8 of [15]), which shows a reasonable agreement. The corresponding calculated hydride volume fractions 
for the three hydrogen concentrations are: C — [0.0055,0.0115,0.0176]. Since the nucleation energy density g„ is the 
driving force for hydride formation, we have listed their calculated values in table [TVl For example, at C = 100 wppm, 
CT^^ = 150 MPa, gn = 554 MJm"'^, and so on. 



X 10" 



0.8 ■ 
0.7 ■ 
0.6 ■ 
0.5 ■ 
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0.3 ■ 
0.2 ■ 
0.1 ■ 





1000 2000 3000 4000 5000 6000 7000 8000 

Time (s) 



700 



Temperature (K) vs. time (s) 



300 




2000 4000 6000 8000 



1000 2000 3000 4000 5000 6000 7000 8000 

f(s) 



FIG. 4: Results of calculations of the hydride volume fraction {() and the orientation parameter (Or) as a function of time at 
a constant applied stress of 150 MPa in Zr-2.5Nb alloy. The sample contained 100 wppm hydrogen. The inset figure displays 
the cooling history during reorientation. 
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TABLE III: Calculated orientation parameter 6r for hydride formation at various total hydrogen contents C and applied stresses 
a^^ a.t T = 373 K. 



Or (-) 




C (wppm) 




cr^^ (MPa) 


100 


200 


300 


100 


0.0765 


0.2381 


0.4067 


150 


0.5052 


0.8071 


0.8461 


200 


0.9153 


0.9772 


0.9797 


250 


0.9912 


0.9980 


0.9982 


300 


0.9991 


0.9998 


0.9998 



TABLE IV; Calculated nucleation energy densities gn for hydride formation at various total hydrogen contents C and apphed 
stresses a^itp at r = 373 K. 



g„ (MJm-^) 




C (wppm) 




a^^ (MPa) 


100 


200 


300 


100 


552 


540 


537 


150 


554 


541 


540 


200 


555 


543 


540 


250 


556 


543 


541 


300 


558 


545 


542 



IV. DISCUSSION 



The behaviour seen in Fig. |4] may be interpreted as an expression of supercoohng effect, whereupon a metastable 
configuration (local minimum) appears before the system settles dovirn to a more stable configuration. This effect, 
despite its importance, to our knowledge, has not been studied experimentally in a quantitative fashion for the Zr-H 
system. In our study as in the earlier investigation , we assumed that stress orientation of hydrides occur primarily 
during the nucleation stage of precipitation. The kinetics of nucleation is characterized by the nucleation rate. The 
steady-state nucleation rate J gives the number of (hydride) nuclei forming per unit volume and time, expressed as 
[S^ . J — ZvcN exp{—G* /kBT), where Z is the Zeldovich factor (« 0.1), i^c the frequency factor, i.e. the rate atoms 
add to the critical nucleus, N the number of sites available for nucleation, and G* is given by Eq. Note that 

G* (X 1/5^, so a slight increase in gn gives a significant rise to the nucleation rate. 

Nucleation experiments similar to that of Tanaka et al. Q made on the Fe-N system, are valuable to determine the 
hydride nucleation rate as well as the degree of hydride orienting as a function of hydrogen concentration and applied 
stress. Such experiments would allow to identify or verify some of the model parameters used in the computations 
here. 

Finally, an issue worth addressing is the role of late-stage growth or coarsening on hydride stress orienting. The 
late-stage growth follows the Lifshitz-Slyozov law [s^, expressed as 

a3~ag = ai3(<-to), (22) 

where a is the mean particle (precipitate) radius at time t, oq the mean particle radius at time to when coarsening 
commences, and is the rate constant for diffusion-limited coarsening. Such a behaviour has been observed for 
growth of hydride platelets in Zr-2.5Nb alloy [3^. It can be shown that the maximum rate of particle size increase 
occurs at a^ax = 2a, and particles with sizes a < a will disappear with relatively high shrinkage rates. 

Following the arguments of Puis [T^, based on a model by Sauthoff "sTl, we consider a tubular specimen of Zr-2.5Nb 
containing hydrides, which at time to after nucleation half of the hydride platelets are radially oriented and the rest 
are circumferentially oriented across the tube wall. If the major radii of these platelets are denoted by and ag, 
respectively, Sauthoff obtained a relationship between these radii using the Gibbs-Thomson-Freundlish relation, 
namely 

ar = — ^ (1 - g), (23) 



where e is the eccentricity (ratio of minor to major radii), 7 = 7/4- 2e7e and g — a^/ ag. By setting g = 1/2, Eq. (|23|) 
gives a minimum mean radius for which even the largest circumferentially oriented hydrides will no longer grow. 
Using the data in the foregoing section with a^p^p = 150 MPa and e = 0.23, we obtain = 1.1 x 10~^ m. 
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In order to attain a complete orienting in radial direction, all the circumferential hydrides must disappear in addition 
to half of the radial ones. Hence, there would be a factor of four decrease in the total number of hydride platelets. 
Since during coarsening the total precipitate volume remains practically constant, this decrease corresponds to a factor 
of (4)1/3 ~ 1.6 increase in the mean platelet radius. 

We may calculate, for the case of Zr-2.5Nb discussed in the foregoing section, the time needed to increase the 
minimum platelet radius by a factor of 1.6. Boyd and Nicholson [s^l have calculated the coarsening rate constant ao 
for disc-shaped precipitates, viz. 

Using Ceq — Cf^iT) and Da in tableU and other data listed previously, Eq. (|24)) yields ao — 2.715 x 10~^^ m^s"^ at 
T = 523 K. Thus to increase the minimum platelet radius of dr = 1.1 x 10^^ m by a factor of 1.6 would take about 15 
s at T = 523 K, according to Eq. ([22]) . Therefore stress orienting of hydride platelets would be possible in Zr-2.5Nb 
alloy during coarsening. 



V. CONCLUSIONS 



In this paper, we have presented a generic model for calculation of stress-induced orientation of disc-shaped precipi- 
tates in alloys. The model is applicable to situations where the stress orienting occurs primarily during the nucleation 
stage of the precipitation process. We have identified parameters appearing in the model for the case of hydride 
precipitates in a zirconium alloy. The model can be used to compute the degree of orienting as a function of stress, 
temperature and time. It also can be extended to account for the effect of supercooling during phase transition. 
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APPENDIX A: DIFFUSION MODEL FOR PRECIPITATION 



An idealized configuration for second-phase precipitate metal composite is considered. The hydrides are a regular 
array of parallel cylinders embedded in a metal matrix. We posit that solute precipitation is diffusion-controlled and 
diffusion occurs only in the radial direction (cylindrical symmetry). Furthermore, there is no imposed external force 
on the solute atoms. The solute concentration C(r, t) is a function of space r and time t obeying the diffusion equation 
of the form 



dC Id/ dC\ 
= (rD I 



dt r dr V dr J 

This equation is solved subject to the following initial and boundary conditions 

C(r,t) 
C{r,t) 
dC{r, t) 
dr 



Co 


at 


t 


^ 0, 


Cs 


at 


r 







at 


r 


= rr. 



(Al) 



(A2) 
(A3) 

(A4) 



where R is the radius of the cylinder and Tc that of the cell enclosing the cylinder assumed to be impenetrable. 
The impenetrability (zero flux) condition is equivalent to the requirement that C (r, t) shall have the symmetry of a 
regularly spaced lattice of (Trr^)"^ identical cylinders per unit area on which precipitation occurs 17[. The solution 
for this boundary value problem [13, HI] is expressed as 



Cir,t)~Co 
Cs — Co 



n=0 



Jo{KR)\^ 
(A„r, 



Cn{r), 



(A5) 



where 



Cn{r) = Jo{Xnr)Yo{XnR} - lo(A„r) Jo(A„i?), (A6) 



10 



and the eigenvalues A„ are the roots of the transcendental equation 

Jo{XnR)YiiKrc) - Yo{XnR)JiiKrc) = 0. 



(A7) 



Here Jn{x) and Yn{x) are the Bessel functions of the first and second kind, respectively, of order n. 

Let us calculate the fraction of precipitates w of the excess solute as it evolves with time. This is related to the 
total number of solute atoms remaining between R and r-c, expressed as 



Mt^2n r[C{r,t)-Co\dr. 



(A8) 



Denoting the quantity of solute in the metal matrix after infinite time by Moo, then the precipitated fraction of 
particles is w = MtjM^. Combining Eqs. (jA5|) and (jA8[) . we write 



w — \ 



R' 



1 - 



MKR)\^ 



-DXtt 



(A9) 



which is essentially the result obtained by Ham [T3|. Considering now Eq. (IA9|) . we note that each eigenmode decays 
exponentially with time with a different decay time r„ = r^/a^D. Since the large-n eigenmodes decay more rapidly 
with time, only the most slowly decaying eigenmode remains in the long-time limit, namely 



w = 1 - 4 



alil-Ryrl)Y'[{ 



'J o{aoR/rc) 
Ji(ao) 



- 1 



-t/ro 



(AlO) 



This relation can be compared with Ham's approximate expression for namely w ^ 1 — exp(— i/ro), with tq — 
rl/aQD and ofg = 2[ln(rc/i?) — 3/5]^^. Evaluations of the precipitate fraction w as a function oiDt/rl made according 
to Eqs. l|A10p . (|A9[) with six term in the sum, and Ham's approximate expression show that Eqs. (jAlOp and (|A9p 
yield almost identical results, whereas Ham's approximate relation gives a slightly lower w for Dt/r^ < 1, see Fig. \5\ 



1 



r /R = 5 

c 




10" 



10" 



Exact, n=0 

Exact, n=5 

Ham's approx. 



10" 10' 



Dt/r 



FIG. 5: Precipitate fraction w for diffusion to a cylindrical sink of radius R surrounded by an impenetrable concentric cylinder 
of radius Vc for rc/R = 5. Here "Exact" n — corresponds to Eq. (|A10|l , and n = 5 to Eq. (|A9[I with the first 6 terms in the 
sum; which are compared with Ham's approximate solution. 
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